Anomalous Diffusion of particles with inertia in 
external potentials 

S. Eule, R. Friedrich and F. Jenko 
February 1, 2008 

Abstract 

Recently a new type of Kramers-Fokker-Planck Equation has been pro- 
posed [R. Friedrich et al. Phys. Rev. Lett. 96, 230601 (2006)] describing 
anomalous diffusion in external potentials. In the present paper the ex- 
plicit cases of a harmonic potential and a velocity-dependend damping are 
incorporated. Exact relations for moments for these cases are presented 
and the asymptotic behaviour for long times is discussed. Interestingly 
the bounding potential and the additional damping by itself lead to a sub- 
diffussive behaviour, while acting together the particle becomes localized 
for long times. 

1 Introduction 

In recent years the concept of Continous Time Random Walks (CTRWs) [T] 
has been extensively used to model systems that exhibit the phenomenon which 
has been termed "Anomalous Diffusion". Various mechanisms in physical sys- 
tems are knovifn leading to sub- and superdifFusion respectively. CTRWs are 
thus suitable to model a variety of physical, chemical, biological and even med- 
ical systems, like charge transport in disordered systems 2 , protein dynamics 
[3], transport in low-dimensional chaotic systems 0], anomalous transport in 
plasmasfSl El E] and the spread of a pandemic [5]. For a review see e.g. [S]. 
CTRWs in a natural way lead to a description within the framework of frac- 
tional kinetic equations and in particular Fractional-Fokker-Planck Equations 
(FFPEs). 

In a recent publication [TU], [TT] the concept of CTRWs has been extended to 
particles with finite inertia under the influence of a stochastic force. In the 
context of this model the particle is subjected to a series of random kicks, such 
that the particle's motion is deterministic most of the time. This leads to a gen- 
eralized equation of Kramers-Fokker-Planck (KFP) type which is characterized 
by a collision operator that is nonlocal in space as well as in time. For a special 
choice of the time-evolution kernel a Fractional Kramers-Fokker-Planck Equa- 
tion has been obtained. Fractional kinetic equations in general and fractional 
diffusion equations, in particular, have a wide range of applications. We should 
mention that over the past years several generalized KFP-equations [T^], [T5] . 
p3] have been proposed (for a review including discussions of various applica- 
tions see e.g., Ref. [H]). In a recent paper [T^ the connection of these equations 
is established unsing an approach based on the underlying Langevin equations. 
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In this paper we want to extend the discussions of the force-free case started in 
[TT] and discuss the effects of an inclusion of a Hnear force and an additional 
damping between the kicks. This case should be of importance for various 
practical applications in the context of anomalous diffusion in random and ac- 
cordingly complex environments. 

For results concerning anomalous transport of particles without inertia in exter- 
nal fields the reader is refered to [T^ . Recent results on CTRWs in an oscillating 
external field are given in [T^. 

This paper is organized as follows. First we present the model under considera- 
tion and incorporate the effects of the harmonic potential and the damping. As 
we will see, the resulting equation will be rather complicated. We thus concen- 
trate on the equations for the moments and discuss in this context some limiting 
cases. We will present analytical and numerical results for the time-evolution 
of the second-order moments. 



2 The model 



Some 70 years ago Kramers investigated the motion of Brownian particles in 
a fluid. In a seminal paper [19j he considered the joint-probability-distribution 
of position and velocity /(x, u,t) for which he could derive the well-known 
equation: 

'd_ 

di 



-u + Vu- A(x, u) 



/(x,u,i) = £Fp/(x,u,t) 



(1) 



where A(x, u) is the acceleration due to an external potential and Cpp is the 
Fokker-Planck collision operator [501 , [5T] 



£fp/ = 7V„- (u/)+i^A„/. 



(2) 



Note that this notation indicates that the divergences act on the product of the 
following function and the probability distribution. 

In a new type of the Kramers-Fokker-Planck Equation has been proposed. 
In the context of this model the motion of the particle is deterministic while 
from time to time it is subjected to a kick in a random direction with random 
velocity. As the time r between the kicks as well as the new velocity u is given 
by some probability distribution, a CTRW-modcl in the velocity-space has been 
obtained. As a main result from the following master equation 

d 

— + V;j; • U + V„ • A(X, U 

ot 



/(x,u,t) = 
^ dt'<S>{t-t') y"du'F(u;u')-P*'*7(x,u',i') 
- [ di'$(i-t')^*^*'/(x,u,t') 

"'0 



(3) 



the generalized Kramers equation 



d 



--+V^-U + Vn- A(X, U) 

ot 



/(x,u,t) = 

^Cfp f dt''S>{t~t')V'''fi^,u,t'), 
Jo 



(4) 
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has been derived, where F{u; u')du' denotes the probabihty that the particle's 
velocity will end up in the velocity space element du about u and ^{t ~ t') is 
some memory kernel. The deterministic propagation between the random kicks 
is governed by the Frobenius-Perron operator 

7'*'*7(x, u, t') = e-(*-*')[^-"+^-^('''")l /(x, u, t') . (5) 

Eq. ([4]) generalizes Eq. ([T]) in two aspects. First it incorporates effects which 
are nonlocal in time and second it contains retardation effects due to P*'* which 
renders the equation nonlocal in space as well. For the special case of an 
asymptotic power-law waiting-time distribution ^(r) ~ 

T--(i+'5) one obtains 

/q ^{t — t') ^ oT^t^^ where o^^t is the Riemann-Liouville fractional deriva- 
tive. For an excellent account on CTRWs the reader is referred to [7] . A detailed 
derivation of Eq. Q with an explicit definition of F{u; u') is found in [TU] . 

3 Damped motion in a harmonic potential 

In this section we shall consider the special case, where the external potential 
is given by A(x, u) = — 2r/u — WqX. The corresponding equations of motion for 
an individual particle without the impact of fluctuations read: 

±(t)^u{t), u(i) = A(x,u) = -277U-tJoX. (6) 

This set of equations determines the behavior of an individual particle between 
two succesive kicks. As we have already mentioned the effect of the Perron- 
Frobenius-Operator is to project the position and velocity of the particle back 
to the last kick. For this reason we have to solve Eq. ^ and invert this solution 
to retain the retardation of the probability distribution function (pdf) on the 
right hand side of Eq. (j4|) . The solution reads 

x(i) = e-"* (^xo coa{iu'{t - t' j) + "° sm{oj'{t - i'))) 

u(<) = e""* l^ -V^o-{v^+u;'^)^o _ ^ _ 

Hereby we imposed the initial conditions x(i') = Xq, u(t') = Uq and used the 
abbrevation uj' — \/oJq — rp . Note that for rj > ujo, uj' becomes purely imaginary 
and we have an exponential relaxation of x(t) and u{t). Inverting this equations 
gives a rather complicated expression and solving Eq. with this retardation 
would be a hopeless task. In the remainder of this paper we thus first restrict 
ourselves to the special cases where we have either set Wq = or one of the 
damping constants ry = 0, 7 = respectively. For these special cases rigorous 
results are presented in the next sections. 

4 Damped motion between the kicks 

In this section the limiting case loq <ti t] is studied i.e. we will analyze purely 
damped motion between the kicks. The motion of such a particle in phase-space 
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is visualized in figure [TJ The motion of an individual particle between the kicks 
is given by 

x(i) = u{t) , u{t) = A(u) = -277U , (8) 

with the solution 

x(t) = xo + Uo^(l-e-2')(t-*'; 

u{t) = Uoe-2')(t-*') . (9) 

As already mentioned we have to invert Eq. ^ to get the retardation of the pdf 
in Eq. (HI, i.e. we have to regard the effect of the Perron- Frobenius-operator 

xo = x(t)-u(t)JL(i_e2''(-*')) 

Uo = u(0e2''(*-*') . (10) 

To calculate the temporal behavior of the moments we use Eq. Here we 

have to be careful. On the lefthand side of (jl]) we have to average over the 
values at t while on the righthand side we average over retarded values at t' . 
Consequently we get retardation effects in the moment equations. Assuming 
the system to be infinite and spatially homogeneous we obtain for the moments 
of lowest order 



d_ 

di 



q{t) = Aq{t) + f $(t - t') (Bq(t') + I) dt' . 
Jo 



(11) 



Here we have introduced 

q = f (ux) I , I - I I , A = I -2?7 1 I and (12) 
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B = -7| e-2'?(t-t') e-2'?(*-*') (^1 _ e-2'j(t-*') j |. (13) 

It is straightforward to generalize these equations for moments of any order. 
Due to the occurence of the convolution integrals it is appropriate to switch to 
the Laplace domain, where this set of equations can be solved. 

<"^)(^)- A + 4,tT,L(2,.A) - (^^) 

= A + 2, + 7C>(2, + A) ^''^ 

(x2)(A) = ^(ux)(A). (16) 

Hereby we have used the Laplace shifting theorem. It is apparent how the ad- 
ditional damping affects the time evolution kernel. Obviously the asymptotic 
behaviour for large times , i.e. the small A-behaviour, is given by the Laplacein- 
version of $(A) for (u^)(A) and (ux)(A) and by for (x^)(A) respectively. 
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To obtain concrete results we consider the important case of a fractional evolu- 
tion kernel, i.e. $(A) = A^^'^. We get 



- , .l^M-. A^^ (17) 



2?7 + 7(2?/) 

which corresponds to the asymptotic behaviour 

K 



r((5)(2ry + 7(2?7)i 



V^^AtJ 4r/2 + 7(2r;)2-« r(^)(277 + 7(277)1-*^ ' ^ ^ 

For the mean-square-displacement we obtain 



, .w. ^ 27((M!:!±M^ ^ (20) 



The additional damping thus shifts the subballistic-superdiffusive behaviour to 
the subdiffusive regime. From Eqs. (fTijl - lfTHI it is clear that the magnitude of 
the damping constant is crucial for the rate of convergence to the asymptotics. 
In the classical limit, i.e. 5 — 1 standard diffusive behaviour {x^){t) ^ t is 
obtained. 

To complete our discussions we mention that for 77 ^ we get the results 
described in [TT]. For 7 ^ we obtain 

We see that even the additional damping alone leads to subdiffusive behaviour. 
This is in contrast to the free flight case, where even if 7 7^ superdiffusive- 
suballistic behaviour is obtained. 



and 



5 Response to an external harmonic potential 

Let us now examine the motion of the particles exposed to an external linear 
restoring force. In other words we will study the effects of an external harmonic 
potential acting on the particles between the kicks. As mentioned before dis- 
cussing this problem with both dampings acting on the particle leads to bulky 
equations. For this reason we consider first the situation when 77 <C Wg. Eq. ([7]) 
then reads 

x(i) = xoCos(cj(<-<')) + — sin(Lj(t-t')) 

UJ 

u{t) = -ujxosin{uj{t-t')) + uocos{uj{t-t'j), (24) 
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Figure 1: Visualization of damped motion between the kicks in phase-space 



where we have imposed the initial conditions x(t') = Xq and u(t') = Uq. In this 
context a word of caution shall be appropriate. One has to be careful in setting 

the initial conditions. 

We can now repeat the procedure of the previous chapter for the calculation of 
the moments. The non-locality in space provides for a mixing of the space and 
velocity coordinates. 
We obtain 

d_ 

di 



q{t) = Aq(i) + f ^{t- t') (Bq(t') + I) dt' . 
Jo 



(25) 



Here we have introduced 

^2 




and (26) 



/ \ 

-7 sm{2oJo{t - t')) cos(2wo(i - i')) ^ ^H'^^oit - t')) 

^ 2(jlsm'^{(Jo{t-t')) -2(Josm{2(Jo{t-t')) cos'^{(Jo{t - 1')) j 

(27) 

Switching to Laplace space we obtain a linear algebraic system 



B 



Aq(A) =C(A)q(A)+I (28) 



where 



C= ^iA-B)-u;o -UA + B) -^iA-B) + l 

V 2^(A + B)-7c.o'*(A) ^A-B)-2u;l -X{A + B)+^^X) 

(29) 

and I = ^0,0, 2X-^^y Hereby we have used the short-hand notation A = 
$(A - 2iwo) and B = $(A + 2iujo)- Note that the relationship ^{z*) = ^*{z) 
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Figure 2: Visualization of the motion in a harmonic potential 



holds and thus A* = B, where the (•)* means complex conjugation. It follows 
in this case that C has only real-valued entries. 

In the Laplace-domain we obtain a linear system of equations which is solvable 
but we have not yet been able to transform these expressions back or to extract 
useful limiting cases. One should notice that for this system the knowledge of 
just the long-time behaviour is no longer sufficient. In particular the oscillations 
of the moments for small t are ignored. Employing the fractional time kernel 
we get for 6=1 the ordinary Kramers equation. It is interesting to note how 
in this case the additional couplings in the dynamical system ([25|) disappear. 

Let us now concentrate on the general case for 7 = 0. The equations for the 
moments up to second order are then no longer different from a retardation free 
version of Eq. ^ (i.e. the following results for the moments stated here also 
hold for the Kramers equation proposed by Barkai and Silbey [12] )■ This allows 
us to obtain concrete results on the behaviour of the moments. 
It is important to first consider the physical meaning of 7 in Eq.([2]). In the clas- 
sical Kramers equation the damping operates between two successive random 
events. In the generalized Kramers equation Eq.Q the damping acts during 
the kicks. A positive 7 affects the probability F{u, u') in a way that it is more 
probable to have smaller absolute value of u' after a kick. This effect is visual- 
ized in figure [3) Without this damping the energy of the particle would tend to 
infinity for the force-free case. If we have an additional damping between the 
kicks in the generalized Kramers eq. though we can set 7 = without loss of 
physical significance. Thus the limiting case 7 ^ is not solely of academic 
interest. The phase-portrait of such a motion is shown in figured 

The set of equations for the second order moments reads in Laplace space 



(u^)(A) 



A(2?7 + A)(A2 +4?7A + 4w2) > 



(30) 



(ux)(A) 



(31) 



(2r7-f A)(A2+4?7A + 4l^2) ' 
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Figure 3: 



The effect of 7 in Eq. ^ visualized for free motion between the kicks 
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Figure 4: The combined effect of the harmonic potential and damping in phase- 
space 
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Figure 5: t] — 0, loq — A 



(x^)(A)^ .7.^^45^%™. (32) 
i^')W=i'''1fj''\'''fi^ (33) 
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- (2, + A)(A2+4,A + 4c.§) ' (34) 

> - (2r; + A)(A2+4ryA + 4^2) • (35) 

The inverse Laplace- Transforms of these expressions can now be exphcitly cal- 
culated. 

For the mean-square-displacement (MSD) (x^)(t) we obtain 

(x2)(i) = (2-'5-ie-2t('?+v^)Xt<5 (^2e2*v^(r(5, -277<) - T{S)){r]t)-^ + e^*^^ (36) 

(i(V^ - v))-Hm - r{S, {2tV^ - rj))) +i~t{V^ - ?7))-'(r(<5) - r(<5, i-2tV^ ~ 77)))) ) 

Hereby ili is given by a) = 77^— cjg and r(Q!, /3) is the incomplete Gamma-function. 
(The occurence of these incomplete Gamma-functions, which are integral expres- 
sions, is due to the multiple occurence of convolution-integrals when backtrans- 
forming the moment-equations) 

We do not want to state the other complete expressions here but rather discuss 
some qualitative features of the results. Let us in the following focus our atten- 
tion on the mean-square-displacement. We present numerical results for some 
interesting limiting cases and discuss these results. The following results were 
obtained for S = 0.5 and K = 1. 

In figure ([5]) we see the MSD oscillating around the asymptotics ^ . The 
bounding potential thus leads to subdiffusive behaviour. 

The next special case is ujq — 0, i.e. no external potential is present. In figure 
El we see the from ([^0]) expected subdiffusive behaviour. 

Let us next consider the case when the motion between the kicks is over- 
damped, i.e. w'^ < 0. We see from figure ([7]) an exponential decay of the MSD. 
The particle localizes at the attracting center. 
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Figure 6: wq = 0, r/ = 2 
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Figure 7: Overdamped motion between the kicks, r/ = 4, = 2 
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Figure 8: Damped oscillations between the kicks, 77 = 1, cjq = 10 



Finally we consider the case w'^ > 0. This corresponds in the classical picture 
to damped oscillating motion between the kicks. 

The behaviour of the MSD is apparently similar to the behaviour of the clas- 
sical position x{t) of the Harmonic Oscillator. It is quite interesting however 
to see that for 77 = as well as for ujo = the MSD shows subdiffusive be- 
haviour, while when both effects are active the MSD exponentially decays and 
so localizes. In this context we want to mention that the long-time behaviour 
of the second-order moments for the Kramers equation proposed by Barkai et 
al. is calcuable even for 7 7^ and leads to loacalization even if 77 = 0. As the 
Kramers equation of Metzler et al. describes subdiffusion, localization is here 
obtained too. 

Summarizing the bounding force as well as the additional damping turns the 
superdiffusive-ballistic behaviour of the force-free case into the subdiffusive 
regime. If both effects are present the particle localizes eventually at the center 
of attraction. 

6 Conclusions 

We have investigated the generalized Kramers equation proposed in [TU] for 
a particle moving in a a harmonic potential under the influence of damping. 
We presented the appropriate Kramers equation. For a fractional time evolu- 
tion kernel we discussed limiting cases for the behaviour of the moments, i.e. 
ajQ = 0, t ^ 00 and the general case for 7 = 0. For ujq = we obtained a 
transition from the superdiffusive-ballistic to the subdiffusive regime. The rate 
of convergence depends on the value of the additional damping. 
For 7 = we could obtain a closed expression for the mean square displace- 
ment. We presented numerical results for the mean square displacement and 
discussed qualitative features. In this case the harmonic potential as well as the 
additional damping shift the ballistic-superdiffusive behaviour, present in the 
force-free case, to the subdiffusive regime. For ujq = Q this also holds true for 
the general case. Interestingly only both effects together lead to a localisation 
of the particle. 



11 



References 

[I] E.W. Montroll and G. Weiss, J. Math. Phys. 6, 167 (1965) 
[2] H. Schcr and E. Montroll, Phys. Rev. B 12, 2455 (1975) 

[3] W. G. Glockle, and T. F. Nonnenmacher, Macromolcculcs 68, 46 (1995) 
[4] G. M. Zaslavsky, M. Edelman and B. A. Niyasov, Chaos 7, 159 (1997) 

[5] D. dcl-Castillo-Ncgrete, B. A. Carreras and V. E. Lynch, Phys. Rev. Lett. 

94, 065003 (2005) 

[6] R. Sanchez, B. Ph. van Milligen B. A. Carreras, Physics of Plasmas 12, 
056105 (2005) 

[7] R. Balescu, Aspects of Anomalous Transport in Plasmas (loP Publishing, 
Bristol, 2005) 

[8] D. Brockmann, L. Hufnagel, and T. Geisel, Nature 439, 462-465 (2006) 

[9] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000), R. Metzler and J. 
Klafter, J. Phys. A:Math. Gen. 37, R161 (2004) 

[10] R. Friedrich, F. Jenko, A. Baule and S. Eule Phys. Rev. Lett. 96, 230601 

(2006) 

[II] R. Friedrich, F. Jenko, A. Baule, and S. Eule, Phys. Rev. E 74, 041103 
(2006) 

[12] E. Barkai and R. Silbey, J. Chcm. Phys. B 104, 3866 (2000) 
[13] R. Metzler and J. Klafter, J. Chcm. Phys. B 104, 3851 (2000) 
[14] R. Metzler and I. M. Sokolov, Europhys. Lett., 58, (2002) 

[15] W. T. Coffey, Y. P. Kalmykov and J. T. Waldron, The Langevin Equation 

(World Scientific, Singapore, 2004) 

[16] R. Friedrich, S. Eule and F. Jenko, Langevin Approach to Fractional Dif- 
fusion Equations including Inertial Effects, to be published 

[17] R. Metzler, J. Klafter and I. M. Sokolov, Phys. Rev. E 58, 1621 (1998) 

[18] I. M. Sokolov and J. Klafter, Chaos, Solitons and Fractals 34 (2007) 

[19] H.A. Kramers, Physica 7, 284 (1940) 

[20] H. Risken, The Fokker-Planck Equation (Springer, Berlin, 1989) 

[21] N.G. van Kampen, Stochastic Processes in Physics and Chemistry (North- 
Holland, Amsterdam, 1981) 



12 



